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A dynamical model for the collapse of the wave function 
in a quantum measurement process is proposed by consid- 
ering the interaction of a quantum system (spin- 1/2) with 
a macroscopic quantum apparatus interacting with an envi- 
ronment in a dissipative manner. The dissipative interaction 
leads to decoherence in the superposition states of the ap- 
paratus, making its behaviour classical in the sense that the 
density matrix becomes diagonal with time. Since the appara- 
tus is also interacting with the system, the probabilities of the 
diagonal density matrix are determined by the state vector of 
the system. We consider a Stern-Gerlach type model, where 
a spin- 1/2 particle is in an inhomogeneous magnetic field, 
the whole set up being in contact with a large environment. 
Here we find that the density matrix of the combined system 
and apparatus becomes diagonal and the momentum of the 
particle becomes correlated with a spin operator, selected by 
the choice of the system-apparatus interaction. This allows 
for a measurement of spin via a momentum measurement on 
the particle with associated probabilities in accordance with 
quantum principles. 

PACS No(s): 03.65.Bz 



I. INTRODUCTION 

After about eighty years of remarkable success of quan- 
tum mechanics, the understanding of its measurement 
aspects remains poor and debatable. The root of the dif- 
ficulty is that while the state vector of a quantum system 
has a deterministic evolution, the measurable properties 
of the system can only be predicted probabilistically. The 
linearity of quantum equations allows for solutions which 
are a linear superposition of some basis states, each of 
which may correspond to one of the different values of a 
dynamical variable. So the definitive outcome that is ob- 
tained in a given measurement of this variable can only 
be accounted for by the notion of the 'collapse' of the 
wavefunction to one of the basis states contained in the 
wavefunction. Though the question as to which state the 
wavefunction collapses to can be answered probabilisti- 
cally, quantum mechanics contains no mechanism for the 
collapse. Bohr |l[ and von Neumann [|j postulated that 
this process occurs when the quantum system comes into 
contact with an apparatus which must be described clas- 
sically. Moreover, since the state vector can be expressed 
as a linear superposition of basis states in any number 



of ways, the interaction between the classical apparatus 
and the system decides which property would be mea- 
sured. According to Bohm j|] this implies that the clas- 
sical properties as we observe them are contained only as 
potentialities in the state vector. Many authors have ex- 
pressed dissatisfaction with this dichotomy between the 
quantum world and the classical world. Are there really 
two kinds of systems? If so, is there a way to describe the 
combined dynamics of quantum and classical systems to 
exhibit the so called 'collapse'? No satisfactory answer 
to these issues has emerged even after intensive efforts by 
a large number of workers Q-^]. 

An interesting line of investigation to resolve this issue 
was initiated by Zeh |^] , who observed that the measure- 
ment apparatus, being always a macroscopic object, has 
closely spaced energy levels which make it very suscep- 
tible to the influence of the environment. The environ- 
ment consists of a large number of degrees of freedom 
and its interaction with the apparatus causes decoher- 
ence in the quantum evolution of the latter. This deco- 
herence is quite a general feature whenever one considers 
the interaction between a small quantum system and a 
large one and monitors the density matrix of the small 
system only . Here the off-diagonal elements of this 
reduced density matrix decay in time. However such de- 
cays seem dissipative only for small times as they are 
actually arising due to the superposition of a large num- 
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the condition that ujj are closely spaced, give rise to ap- 
parent decays for time t T, where T is a recurrence 
time which can be astronomically large, even under mild 
conditions of n ~ 1000 and Alu ~ 10~^. Thus a quan- 
tum system coupled to an environment consisting of a 
large number of degrees of freedom behaves like a clas- 
sical system in the sense that at time scales of interest, 
its density matrix is driven diagonal. Since the diagonal 
density matrix is interpretable in classical terms, several 
authors have examined the effect of environmental inter- 
actions and state reduction in a variety of systems and 
circumstances [^-18|- 

In this paper our goal is to further this line of inquiry 
by examining the process of measurement in the following 
way. We consider a canonical spin-1/2 quantum system 
interacting with an apparatus, which in turn is interact- 
ing with an environment. The basic idea of the scheme is 
to study the following two essential features of the quan- 
tum measurement process. The first is the establishment 
of correlations between states of the system and the appa- 
ratus. The second is the reduction of the density matrix 
of the macroscopic apparatus to a diagonal form dictated 
by the correlations of the apparatus with the system. The 
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second aspect is achieved here by interaction of the appa- 
ratus with the environment. Thus through this scheme 
one is able to study the interaction of a quantum system 
with a macroscopic apparatus within the pure realm of 
quantum mechanics and provide a realization of the early 
ideas of Bohr. 

We consider the set up of a Stern-Gerlach apparatus, 
i.e., a spin-1/2 particle in the presence of an inhomoge- 
neous magnetic field. A measurement of spin is made by 
studying the position or momentum of the particle. Such 
a model has been considered in some detail by Bohm [|. 
Here spin plays the role of the system, and the posi- 
tion/momentum degrees of freedom of the particle that 
of the apparatus. The particle is further coupled to an en- 
vironment via its position, and this coupling is intended 
to drive its translational behaviour classical to perform 
a measurement. The effect of the environment on the 
translational degrees of freedom is taken into account 
via the density matrix equation which incorporates both 
the quantum evolution and the stochastic Fokker-Planck 
type evolution arising due to the environmental interac- 
tion. This equation has been obtained in a variety of ways 
in recent literature ^0 13,1^]. We solve this equation ex- 
actly by including an appropriate system-apparatus in- 
teraction and show that in this case the reduced density 
matrix is driven diagonal at long times in the momentum 
space and has the desired correlation with a component 
of spin. 

The remaining paper is organized as follows. Section 
II contains the calculations and analysis of the density 
matrix for the Stern-Gerlach apparatus, and section III 
summarizes our various observations arising from these 
calculations. 



II. DENSITY MATRIX FOR THE 
STERN-GERLACH APPARATUS 

We consider a Stern-Gerlach type set-up for investi- 
gating the measurement of spin. The Hamiltonian of the 
combined system/apparatus and environment is p9[ | 
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Here x and p denote the position and momentum 
(taken in one dimension for convenience) of the parti- 
cle of mass TO, Acr^ the Hamiltonian of the system, e 
the product of the field gradient and the magnetic mo- 
ment of the particle, H^^ the interaction of the envi- 
ronmental degrees of freedom with x, and denotes 
the Hamiltonian for the environmental degrees of free- 
dom. Since the problem of motion of a spinless particle 
in simple potentials and in interaction with the environ- 
ment has been studied at great length in recent literature 
1^,^-15 1, we draw upon this work to deal directly with a 
reduced density matrix equation for the particle in which 
environmental degrees of freedom have been traced over. 
Though the density matrix equation has been derived in 



a number of ways, the derivation which is in the spirit 
of the present work was given by Caldeira and Leggett 
using Feynman- Vernon path integral approach In 
this method the path integral expression for the density 
matrix is written for the Hamiltonian of (Q), and then 
the degrees of the environment which consists of oscilla- 
tors are integrated out. In the limit of high temperature 
(weak coupling) the expression for the reduced density 
matrix is seen to be a solution of the above mentioned 
density matrix equation. This density matrix equation 
can be thought of as a Markovian limit of a Generalized 
Master equation, and thus its validity lies in the large 
time domain, t ^ tc, where tc is a short relaxation time 
associated with the environment. This limit is needed 
here as the earlier work [|5| shows that the solution then 
evolves to a classical stochastic distribution. The tracing 
over environmental degrees of freedom may be construed 
as if one is not dealing with one measurement process. 
However, we feel that this is an essential characteristic 
of the interaction of a system with a macroscopic object 
and may be viewed in the same spirit as a single measure- 
ment is regarded as an ensemble average for macroscopic 
systems. 

We look at the time evolution of the density matrix 
in the \s,x) representation, where \s) refers to the eigen- 
states of (7z and \x) are the position states. Correspond- 
ing to the four elements of the spin space (tt, 11, Ti, it); 
the equations for the elements of the reduced density ma- 
trix pss'{x,y,t) for our Hamiltonian are: 



dpss'{x,y,t) 
dt 
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d_ d_\ _D 
dx dy ) ih 
ie(xs — ys') i\{s — s 
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Pss'{x,y,t), 
(2) 



where s,s' = + 1 (for f) or -1 (for j). Here 7 is the 
Langevin friction coefficient and D has the usual in- 
terpretation of the diffusion coefficient. In the case of 
a heat bath of harmonic oscillators at temperature T, 
D — 2^mkBT. It is more convenient to work with vari- 
ables: r = [x — y), R = {x + y)/2. Then the spin diagonal 
density matrix pd, and the spin off-diagonal density ma- 
trix Pod obey the equations: 

dpd{R,r,t) ~h d'^pd dpd n 2.fe2 , ^e?" 
ot im oRor or n 

(3) 



where the sign in the last term corresponds to p|| 
and '-' to pii, and 



dpodiR,r,t) _ -hd'^pod 
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where the upper signs in the last two terms correspond to 
p^l and the lower ones to pj^-f. To solve these equations, 
it is convenient to take a partial Fourier transform in the 
variable R: 



piQ,r, t) 



exp{iQR)p{R, r, t)dR. 



(5) 



The equations (^) and ^ simplify to a pair of first-order 
partial differential equations: 

dpd{Q,r,t) /HQ \ dpd Dr^ ier 



dt 



dpod{Q,r, t) 
dt 



dr 



Tt-Q \ dpod 

7f 

m I or 



Dr^ 2e dpod , i>^r 

Such equations, being of first order, can be solved by the 
method of characteristics ||2^. The physical significance 
of the solution can be clearly understood if we choose 
the initial condition to be the following Gaussian wave 
packet of width a and mean momentum p: 



1 



■0(2;, 0) = - — exp(ipa; - I2a^). 



(8) 



The solution for (^, i.e., the spin ofF-diaginal elements 
of the density matrix, for the initial conditions of is 
(see appendix): 



The solution has a factor going as e~^'^^ which drives 
the entire expression to zero with time, independent of 
all other arguments in the density matrix. This means 
that the density matrix is driven diagonal in the spin- 
space. The time scale over which this happens is given 
by 



(11) 



The the solution of (^, i.e, for the spin-diagonal is (see 
appendix for details) 



Pd{Q,r,t) = exp 



+ip{r-rQe - ^(^"^{{r - rgfe - 

(r - rQ)e-^} - i?/(4;iS){r|T + 2rQ 
(r-rQ)(l-e-^) + (r-rQ)2(l-e-2 

Tie/hj{rQT + (r - rQ)(l - e"^)} 
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where rg = hQ/mj. To understand the measurement 
aspect implied by this solution, we consider the solution 
in the momentum representation, i.e.. 



Pd{u,v,t)^ / pd{x,y,t)e'^^-+''yUxdy. 



(13) 



This is obtained by taking a Fourier transform with re- 
spect to the variable r in ( [l2| ) and identifying Q = u — v 
and q = (u + v)/2. This solution is 



PodiQ, r, t) = exp 
exp 



/ T^e^D \ f±2i\t 



4 vri^7 
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eD 
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where t = , and 

rq = Qh/mj ± 2er/7(7,7^ ± 2e/m'-f'^ . 



(9) 
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Pd{Q,q,t) 
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where 

N{t) = {D/2h^-i){\ - e-2-) + (1/^2)6-2-. 



(14) 



(15) 



Now we see that the momentum off-diagonal components 
(Q 7^ 0) vanish with time, reducing pd to the diagonal 
form pO| , p^ -p^ . The time scale over which this happens 
is given by 



777^7^ 



(16) 
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The momentum distribution function can be obtained by 
looking at the diagonal elements of the density matrix 
(p^), i.e., for Q = and q = u : 



[u + pe 



(17) 



This has the classical Ornstein-Uhlenbeck form with the 
spin-dependent drift caused by the field. In the large 
t limit, the momentum distribution is centered around 
e/h'y (—e/h'j) for up (down) spin. Thus we see that the 
measurement of momentum of the particle can determine 
the spin. It is also interesting to consider this solution in 
the space coordinates. The diagonal distribution in mo- 
mentum necessarily implies that the density matrix does 
not reduce to a diagonal form with respect to space coor- 
dinates. The Fourier transform in Q and q of pd{Q,q,t) 
gives the density matrix in position representation: 



Pd{R,r,t) = 2 



M{t) 



exp 



1 

4^' 



D 



(1 



ipe 



-5" 



1 



M(r) 



R 



^(1 

m7 



-2t 



:(1 



-2t 



r) 



ihr 



-2t 



(1- 

where t = and 
M(t) = ct^ + 



iDr 



(1 



2CT^m7 

(1-6-2^)2 



(18) 
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(2r - 3 + Ae- 



-2r\ 



(19) 



As t — > cxD, one can see that the off-diagonal elements 
of the density matrix do not vanish, and the diagonal 
elements give the position distribution function obtained 
by setting r — and R = x: 



|V(x)p=2 



TT / — 1 , _|_ e(l ~ T 

M{t) \ M{t) ^ m-f 



(20) 



The results corresponding to ±e are for the up and down 
spins. The centers of the position distribution function 
shift with time and are clearly different for up and down 
spins. Though the distribution pd{R,r,t) implies a non- 
locality through its dependence on r, we observe that the 
width wi of the distribution in r is considerably smaller 
than the width W2 of the distribution in R. In the large 
T limit, W2/W1 = fi^-m?^^ / D'^ . If pd{R,r,t) is coarse- 
grained over length scales I, such that I is larger than 



the deBroglie wavelength of the particle 7/e and wi, but 
smaller than W2 , we should have a local distribution 
in position space. This means I > Max{'y /e,"f / D), and 
I < DTl^t /vn?^"^ ^ which is surely possible for large enough 

T. 

Thus if the initial wavefunction of the system- 
apparatus is a product of the Gaussian wave-packet of 
(g) and the apparatus state (a| t) + ^1 i)), i-e., 

ij{x,0) = —^-—[a\ ])+b\ i)]exp(-ipx-a;V2a2), 

(21) 

in the model without the environment, the time evolution 
of the density matrix is 

P=I«I'|T)(T |V';(a:,t)^+(y,t) 

+ab*\ T)a \r+{x,t)i>.{y,t) 
+a*h\i){Ur-{x,t),k+{v,t), (22) 

where '>p±{x,t) are the wavefunctions of the particle in 
the potential ±ex. The environment causes the decay of 
the off-diagonal elements and the large time limit of the 
density matrix assumes the form 



PK = |ant)(T Im 



\Pii, 



(23) 



with I and pn being given by ( [l7| ) in the momentum 
representation, and by (^8|) in the coordinate represen- 
tation. This calculation clearly establishes the measure- 
ment of spin via a momentum measurement. The spin 
diagonal density matrix evolves to a diagonal form in 
the momentum space, while the spin off-diagonal density 
matrix goes to zero with time. Further, the probability 
distributions of up and down spins are also given by the 
initial amlitudes a and b according to the quantum pre- 
scription. Figs. 1(a), (b) and (c) show the real part of 
the sum of the density matrices, p = Real{pii + pn), 
given by ( |l7| ) in the momentum representation for the 
diagonal spin elements for three different values of the 
scaled time r — ^t. As r increases, the off-diagonal ele- 
ments clearly decay, leaving a diagonal distribution which 
is centered around two different mean momenta, corre- 
sponding to up and down spins. One can identify these 
mean momenta as those with which the centers of the 
wave packets corresponding to up and down spins move. 

III. SUMMARY OF RESULTS 

In this investigation we have considered through a 
canonical example two essential aspects of the measure- 
ment problem, which are (a) decoherence of the super- 
positions in the apparatus states so as to allow classical 
inference, and (b) the definite correlation of the system 
states with the apparatus states. We find that measure- 
ment is achieved if the apparatus is macroscopic enough 
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to be affected by an environment and furthermore, its rel- 
evant degree of freedom has a classical limit in the sense 
of the Correspondence principle. This is to be contrasted 
with the case in which the relevant degree of freedom of 
the apparatus has a discrete spectrum, because in that 
situation the correlation between the system and the ap- 
paratus states is not achieved. To conclude, we have been 
able to provide a scheme of incorporating a concept like 
" classical apparatus" in a purely quantum formalism and 
demonstrate that a suitable quantum apparatus when 
dissipativcly coupled to an appropriate environment does 
perform a measurement. It is in this sense that we justify 
the concept of Bohr and von Neumann that a measure- 
ment requires the interaction of a quantum system with 
a classical system. 
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APPENDIX A: 

Equation (^) is equivalent to the following set of ordi- 
nary differnetial equations : 



dt _ 

1^ =7(^-^q), 



ds 
dpd 
ds 



(Al) 
(A2) 
(A3) 



with tq = hQ/m'j. The invariants of these orbits with 
respect to s are easily found to be 



and 

h = Pd exp 



/i = (r - rQ)e ^* 
or,,. 2rQ 



(A4) 



27 



7 



(A5) 



Clearly, I2 = w{Ii), where w is an arbitrary function. 
This enables us to write 



Pd{Q,r,t) = w(/i)exp 



4r 



27 



±j\rqt 
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w{Ii) is now determined from the initial condition for 
Pd{Q, 1", 0). One can easily see that 



w {{r - rq)e 



= exp 



ip{r - rq)e 
{{r-rqfe~^'<' + 2rq{r-rq)e-'''} 
(r - rq)e~ 



1 

Ah^\ 7 



27 



)±|(r-rQ)e-^* 



(A7) 



Substituting this in _(JA6| ) gives the result ( [l^ ) for 
Pd{Q,r,t). To solve (0) we first make the transforma- 
tion 

Pod = WeMT2iXt). (A8) 

The equation for W{Q,r,t) is now equivalent to the fol- 
lowing set of differential equations : 

dt 



ds 
dr 

ds 
dQ 
ds 
dW 
ds 



1, 



lir - rq), 
±2e/n, 



Ah' 



-W. 



The invariants for this set of equations are 

/i = h/m'yiQ±2et/h) 

/2 = (r - hQ/jm ± 2e/m72)e-^*. 



(A9) 
(AlO) 
(All) 
(A12) 

(A13) 
(A14) 



The third invariant from (All) is obviously a finction of 
Ii and I2, hence. 



W{Q,r,t)^ f{h,l2)exp 
2e 



2 UiT 
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7 7777 
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/(/i,/2) can now be easily determined from the initial 
condition (||) for p(Q,r, 0): 

fih.h) = exp 



0-2 2ei 2 1 ffhQ 2et\ 

— (g±— ) - ^{y—^^) 

(l-e--*) + 7-e-*±^|V7p{(^ 

2etV 
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exp 
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2h J 7' 
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Substituting for fi^h^h) in ( [aT5| ) gives the result of (||) 



FIG. 1. Plot of the sum of the real part of the 
spin-diagonal density matrices in the momentum representa- 
tion, p = Real{p^^ + pii)/a, given by Eq.(ll), versus dimen- 
sionless momenta u — {q + Q /2)a and v = {q — Q /2)a, for (a) 
T = -yt = 0, e/m'j'^ = 0.0, (b) t = 1, e/m'y'^ = 2.0, (c) r = 3, 
e/m-y'-^ = 2.0, with p = 0.2/cr, D/m^'y^ = cr^, m'y/h = 0.5/a\ 
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